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Abstract 



Diversity patterns of tree speeies in a tropieal forest community are approached by a simple lattice 
model and investigated by Monte Carlo simulations using a backtracking method. Our spatially explicit 
neutral model is based on a simple statistical physics process, namely the diffusion of seeds. The model 
has three parameters: the speciation rate, the size of the meta-community in which the studied tree- 

If^ • community is embedded, and the average surviving time of the seeds. By extensive computer simulations 

we aim the reproduction of relevant statistical measures derived from the experimental data of the Barro 
Colorado Island tree census in year 1995. The first two parameters of the model are fixed to known 
values, characteristic of the studied community, thus obtaining a model with only one freely adjustable 
parameter. As a result of this, the average number of species in the considered territory, the relative 

^ 2 '. species abundance distribution, the species-area relationship and the spatial auto-correlation function 

O ' of the individuals in abundant species are simultaneously fitted with only one parameter which is the 

JL^I average surviving time of the seeds. 

1 Introduction 
> 

'nI I The endeavor to explain the observed statistical patterns of species abundance, distribution and diversity 

in ecological communities and to reveal the general mechanisms responsible for these patterns motivated 
f^ I a great number of scientific works in theoretical ecology and statistical physics in the last decades. Initial 

^^ ■ attempts in this sense consisted in fitting statistical distributions to species abundance data of various 

^— V I samples without paying attention to the underlying ecological processes [IJE]- Theories counting the 

pg i presence of some governing forces in the development of the structure and biodiversity of ecological 

communities have been also elaborated. The traditional theory of niche differentiation attributes this 
guiding force to competition [SHS] . According to this approach, each species in a given community is the 
best competitor in its own ecological niche, therefore no two or more species can coexist in the same niche 
for a long time. In the view of the niche-assembly theory, species can live together in a community only 
when they exhibit differences from one another in the use of the available resources. This approach proved 
to be successful in describing the distribution and abundance of species in various environments , but 
evidently, it has some difficulties in explaining the richness of species often observed in communities 
comprising ecologically similar species, such as tropical forests. The neutral theory [8l[9], introduced as 
an alternative to the niche theory, proposed a striking aspect: it propagates the functional equivalence of 
trophically similar species and attributes the key role to randomness as a driving force in the development 
of the observed ecological patterns. 

Since its publication, the neutral theory has received much attention. This new approach provided the 
first models to explain species abundance based on fundamental biological processes of species extinction, 
immigration or speciation [9Ullj. The striking similarities with classical statistical physics models mad 
it very attractive for the physics community |12H16j . Despite their amazing simplicity, neutral models 
produced accurate predictions for some taxa Plll2lll5| . However, the assumption of functional equivalence 
of individuals regardless of species in an ecological community has generated serious debates among 
ecologists [171 - 120] . In these days, the neutral approach is more and more received as a null model of 



community dynamics and attempts arc made to bridge it to the niche-assembly theory 

The neutral theory of biodiversity is considered analogous to the neutral theory of molecular evolution 
in population genetics [23]. In fact, the neutral theory in ecology traces back to the theory of island 
biogeography |25| and extends it by taking into account spcciation and the number of coexisting species 
as well. In this simple approach space is modeled implicitly: a local community is considered where 
speciation balances extinction by receiving migrants from a large external source pool (metacommunity) 
and its dynamics is realized at the level of individuals. Besides the functional equivalence of all individuals 
of all the species in a local community, the neutral theory also assumes a constant size of the community 
predicting zero-sum multinomial distribution for the species abundance distribution. 

Application of neutral models to describe the relative species abundance (RSA) distribution of tropi- 
cal rain forest communities turned out to be especially resultful. In HubbcU's model [^[T^ , the ecological 
drift, without any additional mechanisms, is sufficient to produce patterns of species diversity and abun- 
dance in accordance with those observed in nature. In the past decade, a series of neutral models based 
on Hubbell's simple spatially implicit model have been developed [T5 1 [T5i rroi[^BH25] . Recently, models 
with various spatial structures [29l - [3T] and spatially explicit models [Til [321435] have been introduced. 
Such neutral models present the metacommunity in a more realistic way and are able to capture to a 
larger extent the effect of individual birth, death, dispersal and speciation processes. Spatially explicit 
neutral models performed well in the prediction of the observed species-area relationships (SAR) in trop- 
ical tree communities. While the simplicity of spatially implicit (mean-field type) models facilitates their 
analytical treatment, the more complex spatially explicit models exhibit limited analytical tractability. 

In the present work a spatially explicit neutral model for tropical tree diversity patterns is introduced 
and studied. There are a great number of ecological patterns worthy of note but reproduction of all these 
patterns simultaneously exceeds the ability of current modeling approaches. However, the importance of 
fitting multiple patterns at the same time is undisputed as this could better reveal the correct contribution 
of ecological processes to shaping real ecological communities. The motivation of this paper comes thus 
from the endeavor to construct a simple spatially explicit neutral model that, besides the most studied 
RSA and SAR measures, can effectively capture another experimentally verifiable macroecological mea- 
sure, the spatial distribution of individuals of a species, characterized by e.g. the spatial autocorrelation 
function (SAF). 

In the literature there are spatial models that can reproduce the observed species abundance distri- 
bution together with the species-area scaling, however, the applicability of these models to reproduce in 
the same time the observed spatial distribution of individuals of a species generally is not investigated or 
attempts concerning this purpose failed. Our previously developed spatially explicit neutral model j34j 
also shows this limited efficiency. While it is successful in reproducing the measured RSA and SAR in 
tropical tree communities, it produces exponential-type decay for the spatial autocorrelation function 
(SAF) of abundant species instead of the expected power-law type decay. 

Our model is based on the diffusion of seeds [36] , the diffusion assumed to be a simple random walk in 
the two dimensional space. The maximal surviving time of seeds defines the scale of the dispersal kernel, 
which describes the probability for an individual of a species to occupy locations at different distances 
in the ecosystem. For dispersal rule, one can find various assumptions in the literature [I4 ] [32 ] [33 | I36 | [37 ] . 
As already pointed out in [35], the choice of the dispersal distance is crucial in seed diffusion models 
as it strongly infiuences the shape of the obtained species-area curves: on a log- log scale, models with 
nearest-neighbor dispersal result in convex curves, while models with global dispersal produce concave 
curves. 

The approach considered here is based on a three-parameter model, one of these parameters being 
a freely adjustable variable: the average surviving time of seeds (common for all species in the local 
community). The other two parameters are set to realistic values corresponding to the investigated 
system. However, mention must be made of the incertitude of the speciation rate: since its value is 
difficult to measure in practice and in neutral models one can find values for this parameter taken from 



a wide range, the value of the speciation rate set here may also seem speculative. The size of the 
mctacommunity is fixed aceording to the area of an extended region clearly delimited by natural borders 
(e.g. a region delimited by see, ridge of mountains), surrounding the studied region. In real communities, 
the dispersal of seeds can be influenced or restricted by properties of the natural environment itself. 
Delimiting thus the size of the mctacommunity, the model is adaptable to a spatial region where the 
dispersal process of seeds can be described properly by assuming the same dispersal rules. In this manner 
the model essentially provides a simple one-parameter flt for several macroecological measures of the 
investigated tropical tree community: RSA, SAR, species number and SAF. 

The model is applied here to reproduce the ecological patterns obtained for the Barro Colorado Island 
(BCI) tropical tree census [381139] . the oldest and most studied moist forest plot. A number of individual 
trees (corresponding to the number of trees in the studied region of the BCI plot) are positioned on a 
uniform grid, representing the survey area in the model. Computer simulation of the model is realized 
by using the backtracking coalescence technique, described in [JO] and applied in [531lllll32]j a powerful 
simulation method for neutral models in ecology. This method, consisting of tracing lineages of all 
individuals originating in the survey area backward in time, until they speciate or coalesce, proved to 
be efficient for large and complex spatial structures (even for the case of infinite landscapes). In the 
present study, around the survey area, a region corresponding to the area of the whole BCI is involved in 
the calculations, representing the mctacommunity, and a speciation rate for tropical tree communities is 
assumed in the model. Computer simulations of the model focus on the reproduction of the number of 
species, the relative species abundance distribution, the species-area scaling and the spatial distribution 
of individuals of abundant species in the BCI plot, by investigating the effect of the average surviving 
time of the seeds. Since the model has only one freely adjustable parameter, our goal to reproduce 
simultaneously all the mentioned macro-ecological measures seems quite ambitious. In such case one 
should not expect to obtain better fits for all these measures than the fits provided by different models 
focusing on specified ecological measures, considering only one or just a few of them from the whole set 
of relevant characteristics. Fitting only one or two of these measures can definitely result in considerably 
better fits, however, neglecting completely the reproduction of other important properties. In our view 
a successful model should provide a description for all the measurable macroecological patterns, and we 
consider thus that obtaining acceptable qualitative fits, simultaneously valid for the relevant measures, 
is more important than an excellent fit for one or just a few measures. 

2 Experimental data 

As experimental data, the BCI tropical tree census data is used [38 ] l39 j . Based on this, the relevant 
macroecological measures are defined and compared with simulation outcomes. 

The BCI, administrated by the Smithsonian Tropical Research Institute, a biological reserve from 
1923, is located in the Gatun Lake, Panama. The island is covered with moist forest. Since 1980, a 50 ha 
(1000 m X 500 m) region is regularly mapped, in each census all free standing woody stems with a stem 
diameter at least 1 cm at breast height are identified, tagged and recorded. The whole dataset contains 
information on about 240000 stems of more than 300 tree and shrub species, based on this, accurate 
statistics can be made. Census data is publicly available for several years. 

In this work census data from the year 1995 is used. For this year, 243541 individuals of 302 species 
with at least 1 cm stem diameter can be found in the dataset. Analyzing the density of individuals 
in the 50 ha plot, a less populated region can be identified in the western part of it. corresponding to 
a swampland in the studied area (this pattern is observable independently of census year). To ensure 
homogeneity when calculating the ecological characteristics, only a subregion of the entire plot, limited 
to the eastern 500 m x 500 m section, is considered. This 25 ha region contains 112543 trees and shrubs 
of 273 species. 



3 Model 

We investigate a local community in a spatial region of size L x L embedded in a much larger territory 
(e.g. island, continent, infinite land). Our model follows the idea formulated in j36j with the notable 
difference in the determination of the used dispersal kernel: instead of postulating a dispersal kernel 
without arguing its relation to ecological processes, here it is derived from a real physical process based 
on the diffusion of seeds. Following separately the random motion of each seed is a computationally 
demanding task due to the large number of individuals and seeds present in the considered community, 
therefore a statistical description of the spatial distribution of seeds is considered. 

Our basic hypothesis is that from time to time one individual dies and its place is taken by a new 
individual originating from a randomly chosen seed found in its neighborhood. It is assumed that (i) each 
individual emits the same number of seeds at each time-step, which randomly diffuse in the 2D space by 
a simple unbiased random walk; (ii) the seeds survive W time steps, after which they die and become 
unable to grow; (iii) all the seeds (even those belonging to different species) have the same diffusion 
constant and the same surviving period W. 

Individuals of the studied region with area L x L are placed on a uniform grid of size N x N. The 
schematic picture of the studied region is drawn in Fig. [T] The individuals are labeled by mesh coordinates 
(«, j). In the neighborhood of each individual there can be seeds originating from other individuals, and 
we denote by Qij{k,l) the number of seeds in the neighborhood of individual at grid position {i,j) 
originating from the individual at grid position {k,l). 
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Figure 1. Schematic view of the modeled area. The light-coloured square region represents the studied 
area of size L x L embedded in a much larger territory. Individuals are placed to a uniform square grid 
(N X N), labeled by grid coordinates {i,j). The black circles represents individuals in traced "active" 
lineages, color circles stand for lineages finalized by speciation. 



A simple two-dimensional unbiased random walk is assumed for all the seeds. The Qij{k, I) numbers 
are computed for each pair of individuals, according to the basic statistical properties of simple random 
walks. It is assumed that in unit time a seed can make a unit step defined as the grid constant L/N, 
randomly moving to an individual in the neighborhood of its current position in a similar manner as 



the simple unbiased 2D random walker. This defines the time-unit in the simulations. In such case the 
probability P„(x,j/) to find the random walker at position {x,y) after n number of steps when started 
from coordinates (0, 0) writes as: 

P„(x,y) = — expf-^^±^y (1) 

Assuming that the seeds can survive W time steps, the number of seeds originating from (0, 0) that 
will be in position {x, y) can be calculated as: 

Q^,y(0,0)«c/ — expf- — jdz, (2) 

where r = \J x^ + y^ and c is a constant that quantifies the number of seeds emitted by individuals in unit 
time. Expressing the integral term in the above equation by incomplete gamma functions, the number 
of seeds in the neighborhood of the position (i,j) originating from an arbitrary (k^l) position can be 
computed as 



ijK"^! ') 



r(o,-)-r(o,r2) 



(3) 



where r = \J{i — kY + {j — 0^ ^'^'^ the incomplete gamma function r(0,j/) is defined as: 



r(0,y) = y —dt. (4) 

According to the Qij{k,l) quantities one can select the "successful" seed that will generate the new 
individual in an empty place. Details of this dynamical model are given below. 

The model is initialized as follows. Individuals are positioned on a uniform N x N grid and for each 
individual a lineage is considered. The Qij{k,l) quantities are defined for all the (i,j), {k,l) G N x N 
pairs. Initially all the lineages are traced using the coalescence technique [30]. The dynamics of the model 
is defined by the following rules: 

• In each step one individual is randomly selected from all the lineages being traced at that moment 
and its («,j) grid coordinates are determined. The selected individual is the one that was most 
recently born. The origin of this individual has to be defined. 

• With probability q, speciation happens in the lineage. If speciation takes place, tracing of the 
lineage is stopped and the selected individual is considered to be the root of a new species. 

• With probability 1 — g, a parent is defined for the selected individual. The position of the parent 
is chosen according to the dispersal kernel. Coalescence occurs when two lineages share the same 
spatial position. 

The parent selection mechanism of the last rule is done by respecting the weights given by the 
Q{r) = Qij{k,l) values (see Fig. [2]). Based on these values, a seed originating from an individual at 
(fc, /) is chosen, thus, the individual at (fc, I) which emitted the "successful" seed reproduces itself at the 
position (i, j). 

This selection process realizes the simple unbiased two-dimensional diffusion of seeds in the considered 
area, allowing also seeds originating from distant locations to reach and occupy the vacant position. 

The above defined dynamical model can be easily implemented in a computer simulation. Large 
systems, containing the same number of individuals as a known experimental dataset, can be studied in 
reasonable computational time. 
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Figure 2. The number of seeds Q{r) in tiie neigliborliood of a given position originating from a 
location at distance r as a function of the distance calculated for different values of the seed survival 
time, W. r is measured in units of the grid constant L/N. 

4 Computational details 

The model described above is applied here to study the tropical forest community of BCI. In principle, 
for a fixed community with N x N individuals the model has three parameters: the lifetime of the seeds 
W, the size of the surrounding metacommunity and the speciation rate, q. As experimental data, we 
use census data for the eastern (square like) half part of the whole sampled territory (25 ha) of BCI 
for the year 1995. The total number of individuals in this region, K, is 112543. In order to have the 
number of individuals in the model in accordance with the experimental number of individuals, the 
uniform N x N square grid is constructed with A^=335. This grid, representing the survey area from 
BCI is embedded in the center of a 1130x1130 square grid, representing in a proportional manner the 
whole island of Barro Colorado (1560 ha). During simulation, tracing of lineages is realized in this larger 
grid, the position of parents, selected according to the proposed dispersal kernel, can be positions in 
this larger grid. Simulations have been performed by fixing the value of the speciation rate at a value 
comparable to the mutation rate of a neutral gene. In agreement with the accepted order of magnitude 
in the literature [43], q = 1.0 x 10"^ is set as speciation rate and the maximum seed survival period W 
is varied from 50 to 5000 steps in searching for a W parameter reproducing the experimental value of 
the total number of species {S = 273) and the other observed ecological patterns. As already stated in 
the introduction, our model is essentially a one-parameter model, since fixing the speciation rate and the 
size of the metacommunity, we are left with only one freely adjustable fitting parameter. 



5 Quantities of interest 

Apart of the total number of individuals and species in the considered local community there are several 
relevant distributions that should be compared with the census data. For a detailed description of these 
quantities one could consult our previous work |34j . 

For the spatial distribution of species and individuals two relevant quantities are calculated and com- 
pared with experimental results: the Species-Area Relationship (SAR) and the Spatial Autocorrelation 
Function (SAF) of individuals belonging to the same species. 



The SAR is studied by eonsidering larger and larger territories and counting the number of species 
present within these areas. Average species numbers are calculated on subregions with the same size, the 
SAR curve representing these averages as a function of the considered area 

The SAF, denoted by Cs{r) for species s, is constructed in the following way: a uniform mesh (not 
necessarily the same division as for the mesh of the simulation) is considered on the studied area and in 
each {i,j) cell of this mesh, the number of individuals K^ from the considered species is determined. 
The auto-correlation function for a relative coordinate (p, q) is calculated as 



^:] = {(K^l {Ks))iK^,t,+, {Ks)))^^^ (5) 



Here (Ks) denotes the average number of individuals of the considered species on the mesh sites: (Ks) = 
{K^-)ij. The average of the Cp'^q values for all (p, q) mesh coordinates that are inside a ring with radius 

r and width Ar {r < y/p^ + cp < r + Ar), considering a reasonably small Ar value, leads to our definition 

of SAF: 

{cl:l{r<,/^T^<r + Ar}) 
Cs{r) = — ^ (6) 

Another important quantity that should be also compared with the experimental results is the Rel- 
ative Species Abundance (RSA) distribution. This distribution characterizes the frequency of species 
with a given abundance. Different types of representations are generally used to study RSA. The most 
widespread method to present species abundances is due to F. W. Preston [44], who sorted the species of 
a sample into abundance intervals of consecutively doubling lengths ([1,2), [2,4), [4,8), etc.), and plot- 
ted the number of species found within these "octaves". This kind of plotting is motivated by the fact 
that species abundances can vary over several orders of magnitude and usually there are fewer abundant 
species than rare ones. This way the presence of large statistical fluctuations at the tail of the curve can 
be avoided. A second way of representing the species-abundance distribution is plotting the associated 
probability density function p{v), i.e. the probability for finding a species with a given v abundance. 
p{v) can be derived from the Preston plot by dividing the number of species in each interval with the 
length of the interval, and plotting on a double logarithmic plot this quantity versus the mean abundance 
in the given interval. A third way of representing the species abundance distribution can be realized 
by arranging the species in decreasing order by their abundances and plotting the rank of the species 
versus its abundance (rank-abundance plot) again on a double-logarithmic graph. This representation 
is inspired by several abundance studies in sociology and economies, leading to the general Pareto-Zipf 
distribution. 

6 Results 

Fixing the value of g = 1.0 x 10~^, we study the above discussed quantities of interest as a function of 
W, and compare the results with the ones computed from the BCI data. 

In Fig. |3] results for the SAR curves are presented. Simulations are performed with W values ranging 
from 100 up to 3000. With increasing W values, the experimental SAR curve is better and better 
reproduced by our model. 

In Fig. |4]the SAF for the most abundant species obtained with different values of the W parameter 
are presented. For a visual comparison, the SAF for the 1**, 2"*^, 3'"'', 5*'* and 10*'' most abundant species 
in the BCI census, Hybanthus prunifolius (17499 individuals), Faramea occidentalis (12331 individuals), 
Trichilia tuberculata (8684 individuals), Alseis blackiana (4121 individuals) and Capparis frondosa (2078 
individuals) respectively, are also plotted. Simulation curves are averaged results of 100 runs. It should 
be noted that in the simulations, even if the total number of species is the same as in the BCI census, 
the number of individuals of the most abundant species can show considerable fluctuations, this is the 
reason why the simulation curves are vertically shifted relative to the BCI results. 
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Figure 3. Species-area curves. Simulation results (different lines or lines with open symbols) are 
presented for various values of the W seed surviving period with fixed speciation rate [q = 1.0 x 10~^). 
The continuous line with filled circles shows the species-area relation (SAR) for the Barro Colorado 
Island (BCI) census data (census year: 1995). 

From the SAR and S AF results one can learn that the spatial distribution of the species and individuals 
are comparable to the experimental ones for W ^ 1500 steps. The RSA results for W = 1700 are 
presented in Fig. [SJ All three accepted representations of the RSA distribution are plotted and the 
results are compared with those obtained for the BCI census data of year 1995. 

We conclude thus that for W ~ 1500 steps and q = 1.0 x 10~'^ parameters the model generates 
distribution of individuals and species in agreement with the experimental results offered by the BCI 
dataset. Moreover, the total number of species in the system (local community) is also in agreement with 
the number of species in the original BCI dataset. 



7 Discussion 

A spatially explicit neutral model was elaborated to reproduce the observed ecological measures for a 
tropical tree community. Monte Carlo simulations were done on a square-like territory embedded in a 
larger square-like metacommunity, considering a uniform square lattice topology for tree positions and 
using the backward simulation technique. 

In the model two main processes are responsible for the diversity of species: the dispersal of seeds 
capable of growth and the speciation, the parameters governing these two processes are W and q, re- 
spectively. The dispersal of seeds is interpreted here as collection of individual and independent random 
walks on the 2D plane. The W parameter controls the maximum possible dispersal distance (surviving 
time of seeds) and reficcts the importance and validity limits of long-range interactions. The q parameter 
explicitly specifies the rate of speciation within the considered local community. 

The studied area was represented by a uniform grid of 335 x 335 positions of individuals, embedded 
in a larger grid of 1130 x 1130 positions. This larger grid models the metacommunity, which in our case 
is the whole island. The number of trees in the whole island is not known, therefore the size of the larger 
grid was selected proportionally with the size of the BCI. By fixing the speciation rate to a reasonable 
q = 1.0 X 10^^ order of magnitude, we remain with a one parameter {W) model. 

Initially a maximum survival distance of 100 steps was assumed, corresponding to a trail of length 
150 m. For this initial combination of parameters the experimental SAR is not reproduced, there are 
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Figure 4. Autocorrelation functions of tlie most abundant species. Simulation results (lines with 
symbols - averaged curves) are presented for different values of the W seed surviving period with fixed 
spcciation rate (g = f .0 x 10~^). Results obtained for the five most abundant species in the Barro 
Colorado Island (BCI) census data are shown with lines (census year: 1995). 



large discrepancies especially in the limit of large areas. Simulations yield smaller number of species than 
expected, which suggests that the intermixture of species in not effective enough. This means that the 
limit of long-range interactions (governed by W) is farther away than 150 m, which was initially assumed. 
In Fig. |3]the SAR curves with maximum survival distance between 100 and 3000 steps, corresponding 
to trail lengths between 150 m and 4491 m are presented. For maximum survival distances above 1500 
steps, corresponding to a physical distance of a 2247 m, the SAR curves from the simulations exhibit the 
typical power-law like behavior of the experimental SAR curves in the limit of large areas. Increasing 
the maximum survival distance facilitates the intermixing of species, resulting in simulation curves more 
and more close to the experimental one. For W — 1500 steps, comparison of the simulation outcome and 
experimental SAR curves shows 22 % difference of the results. However, one has to be aware of the SAF 
of abundant species as well. For the abundant species of the BCI census area, this function shows power- 
law type decrease for not too large distances. This power-law type decrease was not well-reproduced 
by the earlier spatially extended models. In Fig. H] the variation of the SAF is presented for different 
survival distances in comparison with experimental SAF curves. Analyzing the simulation results for 
values of the survival distance W > 1500 steps, the power-law decrease is reproduced, in agreement with 
the experimental SAF curves. Following the third measure of interest, the RSA, as a function of W, the 
W = 1700 steps case proves to be the most appropriate to reflect the characteristic ecological measures 
of the BCI plot (Figure [SJ. For this case, the average difference in the simulation and experimental 
results is quite high: 53 %. This large difference should be judged by taking into account that wc fit 
simultaneously four relevant ecological measures with only one scalar parameter. 

The results obtained for W — 1700 and q = 1.0 x 10~^ are acceptable thus for the ecological measures 
considered here. Naturally, none of the fits for the relevant ecological measures is perfect, one cannot 
expect this from a simple one-parameter neutral model. Real tropical tree communities arc complex 
systems, and nowadays it is expected that neutral models should be viewed as null models of community 
dynamics, and niche-assembly elements are also needed to obtain a better description of the systems. 
The non-perfect description offered by this simple approach confirms this view. 
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Figure 5. Relative species abundance (RSA). The RSA is plotted in different ways: (a) Preston-plot, 
(b) rank-abundance plot and (c) probability density function. Simulation results are presented for 
W ~ 1700 steps seed surviving period and q = 1.0 x 10~^ speciation rate. Experimental data obtained 
from the Barro Colorado Island (BCI) dataset is also plotted (census year: 1995). For guiding the eye 
in panel (c) a power- law with exponent —1 is also indicated with a dashed line. 

8 Conclusion 

The statistical patterns observed in a tropical tree community were theoretically approached by a spatially 
explicit neutral model based on the diffusion of seeds. The considered model has three parameters, the 
maximum diffusion distance of seeds, the probability of speciation in the system and the size of the 
mctacommunity. Two of these parameters (the size of the metacommunity and the speciation rate) 
were fixed according to the experimental conditions. In principle, we remained thus with a simple 
one-parameter model and this parameter (W) was adjusted in order to reproduce the major ecological 
measures obtained from the BCI tropical tree census. Considering the same number of individuals in 
the model as in the BCI census, we were searching for a proper W value, which would simultaneously 
reproduce the observed equilibrium number of species and the qualitative trends in the RSA, SAR and 
SAF measures. While RSA and SAR are widely studied in the literature and a number of models focus 
on their reproduction, the spatial distribution of individuals of a given species draws less attention. For 
the BCI plot and for not too big distances, the SAF curves exhibit power-law decay. This is true for 
species with various abundance, ranging from a few hundreds of individuals up to the case of the most 
abundant species. The present study considers this measure equally important and besides the study of 
RSA and SAR, the spatial distribution of a given species (usually the most abundant species is chosen 
for better statistics) is also closely investigated. 

Our previous attempt }34j towards a more precise description of tropical tree communities failed in 
the reproduction of the spatial distribution of individual species by assuming constant probability of 
long-range interactions and favoring first neighbor interactions of individuals placed on a uniform grid. 
The present model comprises in a more realistic way the long-range interactions by assuming random 
diffusion of seeds. Also, the simulation method used for studying the model is improved by considering 
the coalescence technique which is faster and makes it possible to take into consideration the role of the 
surrounding metacommunity in determining the species diversity. 

Using a speciation rate of g = 1.0 x 10^^ and assuming a maximum survival distance of about 2500 
m, the model produces results in agreement with the experimentally investigated ecological measures. 
The species-area scaling and species abundance distribution is reproduced in an acceptable manner, the 
autocorrelation function of the most abundant species shows the observed power-law decay in the limit 
of small distances as well. The strength of the model is that with only one freely adjustable parameter 
wc arc able to qualitatively approach a quite complex phenomenon. Obviously, the model neglects many 
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existing ecological processes, however, the results obtained so far suggest that this simple neutral approach 
captures many important aspects of tropical tree communities. 
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